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Resumen 


En este trabajo se presenta un método geoestadístico para realizar estimaciones 
espacio-temporales y se aplica a la estimación de la carga hidráulica del acuífero de 
Querétaro-Obrajuelo en el periodo 1981-2004. Las estimaciones para los años 1993, 
1995 y 1999 se comparan con las de los métodos de krigeado ordinario y krigeado 
multivariado, utilizando el método de validación cruzada, reportándose los valores 
del error medio (EM), error cuadrático medio (ECM) y error cuadrático medio 
estándar (ECME). El promedio del error medio para los tres años elegidos, para 
el caso de krigeado ordinario, resulta ser el menor (-0.23); el valor promedio del 
error cuadrático medio (ECM) es más bajo en el caso del krigeado espacio-temporal 
(224.29), y el valor promedio del error cuadrático medio estándar (ECME) es mejor 
para el método de cokrigeado, ya que sus valores de ECME son cercanos a la 
unidad (0.95). El ECME para el método espacio-tiempo es de 0.8 al considerar todos 
los tiempos; pero en lo particular para los años elegidos, este índice disminuye. 
Los resultados de las varianzas del valor estimado son siempre más pequeños 
utilizando el método espacio-temporal, ya que utiliza más información para la 
estimación; asimismo, es posible hacer las estimaciones en todo el espacio para 
todos los tiempos. Se concluyó que la herramienta es poderosa, porque considera 
toda la información disponible al efectuar las estimaciones. 


Palabras clave: modelo espacio-temporal, geoestadística, estimación carga 
hidráulica. 


Introducción 


Los métodos geoestadísticos de krigeado 
han sido ampliamente utilizados en la 
hidrogeología subterránea. Numerosos inves- 
tigadores usan los métodos geoestadísticos 
para estimar campos de transmisividad 
(Ahmed y De Marsily, 1987) 
hidrostática de un acuífero (Delhomme, 1978; 
Gambolati y Volpi, 1979a y 1979b; Aboufira- 
ssi y Mariño, 1983; Rouhani, 1986; Gambolati 
y Galeati, 1987) o la precipitación (Hevesi et 
al., 1992a y 1992b). 


la carga 


La constante en estos trabajos es que la varia- 
ble se estima en el dominio del espacio; sin 
embargo, sesabe que las variables son dinámicas, 
lo cual indica que tienen una componente 
temporal de importancia. En algunos casos, los 
cambios son lentos (e.g.: transmisividad) y en 
otros son rápidos (e.g.: carga hidráulica). 

De forma tradicional, se ha realizado el 
análisis bajo dos concepciones: el punto de vista 
espacial (ASCE, 199a y 1992b) y el punto de 
vista temporal (Asli y Marcotte, 1995; Larocque 
et al., 1998; Angelini, 1997; Lee y Lee, 2000). No 
obstante, la conceptualización conjunta a la 
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que nos referiremos como espacio-temporal 
ha sido ámbito de estudio en los últimos 15 
años y poco estudiada en el caso particular 
de la hidrogeología subterránea (Christakos, 
2000). 

Rouhani y Hall (1989) son los primeros 
investigadores que plantean determinar 
la carga hidráulica utilizando información 
conjunta en el dominio del espacio y el 
tiempo. Los autores muestran cómo calcular 
dicha estimación usando una modelo del tipo 
lineal: 


Cs: (M,h,)= C, (1,) En C; (h,) (1) 


donde C.(Hh) es la función de covarianza 
espacio-temporal; C (h ), C,(h,) son las funciones 
de covarianza espacial y covarianza temporal, 
respectivamente. El modelo lineal presenta 
problemas, ya que la matriz de covarianza, 
para ciertos casos, es singular (Myers y Journel, 
1990; Rouhani y Myers, 1990; Dimitrakopoulos 
y Lou, 1994). 

El objetivo de este trabajo consiste en 
presentar una metodología para estimar 
la carga hidráulica desde el punto de vista 
espacio-temporal, utilizando un modelo 
espacio-temporal tipo generalizado producto 
suma (De laco et al., 2001; De Cesare ef 
al., 2001). En la metodología se enfatiza la 
determinación de la deriva; el cálculo del 
semivariograma muestral; el modelo espacio- 
temporal y la forma de ajustarlo, mediante un 
proceso visual, y por último, en la validación 
del modelo a través del método cruzado 
(validación cruzada). La metodología fue 
aplicada al acuífero de Querétaro-Obajuelo y 
los resultados se comparan, en términos de 
la validación cruzada y los valores de error, 
con los obtenidos al utilizar sólo el método de 
krigeado y multivariado (Mendoza y Herrera, 
2007). 

La propuesta difiere de la de Rouhani y 
Hall en cuanto al tipo de modelo utilizado: 
producto-suma; pero en la estimación se usa el 
mismo sistema de krigeado (krigeado ordinario 
espacio-temporal). 


Área de estudio 


La zona de estudio comprende una superficie 
de 426 km?, que incluye el acuífero del valle 
de Querétaro, ubicado en la porción sur 
oriental de este estado, y el acuífero del 
valle de Obrajuelo, localizado en el estado 
de Guanajuato y colindante con el acuífero 
anteriormente citado (figura 1). 

En el área se encuentran dos corrientes 
poco caudalosas: los ríos Querétaro y El 
Pueblito. La corriente principal es el río 
Querétaro, 
hasta interceptarse con el río El Pueblito. La 


que corre desde La Cañada 


información geológica y estructural indica 
que los valles de Querétaro y Obrajuelo se 
alojan en una secuencia de fosas tectónicas 
(originada por el tectonismo ocurrido en 
el Terciario), por lo que el acuífero está 
limitado por un sistema de fallas sepultadas 
por sedimentos aluviales, flujos de lava y 
depósitos volcanoclásticos lacustres potentes. 
La disposición natural de las unidades 
litológicas agrupadas en los medios poroso y 
fracturado conforman el sistema que define 
el flujo del agua subterránea en la zona de 
estudio, lo que constituye un acuífero del 
tipo libre, comportándose como confinado 
localmente en algunas zonas por la presencia 
de arcilla (Guysa, 1996). 

Del modelo conceptual se infiere que las 
principales zonas de recarga al acuífero se 
presentan al este, por La Cañada; al noreste, 
por las montañas; al norte, en la zona de Jurica; 
al sur se tiene una entrada en el lineamiento del 
río El Pueblito. 

Se consideran dos sistemas de flujo: uno 
local, formado por zonas con elevaciones 
del nivel freático de poca profundidad en 
la parte central del valle, y un sistema de 
flujo intermedio que sustenta casi el total de 
extracción en el valle, formado también por 
medio poroso y fracturado. 

Por otro lado, en la zona límite entre 
los estados de Querétaro y Guanajuato hay 
abatimientos 


presencia de importantes, 


además de que es imposible definir la 
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Figura 1. Área de estudio. 


presencia de un parteaguas subterráneo o 
alguna frontera física que permita separar 
hidrogeológicamente ambos acuíferos (Simuta, 
2005). La información piezométrica analizada 
en el presente artículo comprende un periodo 
de 24 años (1981-2004), que se detallará en la 
sección de resultados. 


Marco teórico 


En este artículo, la carga hidráulica se 
representará como una función aleatoria 
dependiente del espacio y el tiempo. Así, se 
considera un dominio DeR? en el espacio y 
un dominio T<cR' en el tiempo. Una variable 
aleatoria espacio-temporal (VA) Z(x,f) es una 
variable que puede tomar una serie de valores 


(realizaciones) en cualquier localización en el 
espacio (xe D, en un instante de tiempo te T. La 
VA Z(x,t) quedará totalmente caracterizada al 
conocer su función de distribución, que da la 
probabilidad de la variable Z en la posición x 
en el espacio y en el instante t: 


F(x,t:z)=Prob IZ (x,t)<z) Wz, (x,t)e DxT (2) 


Una función aleatoria espacio-temporal 
(FA) se define como una familia de variables 
aleatorias (Z(x£),/(x,t)e DxT|. Una realización 
de la FA Z(x,f) es una colección de realizaciones 
de las VA. 

Kyriakidis y Journel (1999) han definido tres 
puntos de vista para la modelación espacio- 
temporal. 
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Función aleatoria espacio-tiempo 


Típicamente descompuesta en una compo- 
nente de deriva que modela una variabilidad 
suave de la variable aleatoria Z(x,f) y una 
componente estocástica, la cual denota altas 
frecuencias alrededor de la deriva, ambas, en 
espacio y tiempo (Bilonick, 1983; Bilonick, 1988; 
Bilonick y Nichols, 1983; Bogaert y Christakos, 
1997; De laco et al., 2001). 


Multivariada de la variable aleatoria 


En este marco se conceptualiza la variable 
aleatoria Z(xt) como una colección finita 
de observaciones  Z(x),  correlacionadas 
temporalmente (Solow y  Gorelick, 1986; 
Comegna y Basile, 1994; Papritz y Flúhler, 1994; 
Mendoza y Herrera, 2007). 


Multivariada de serie de tiempo (S/T) 


Modelos de series de tiempo, S/T, son 
conceptualizados como una colección finita 
Z(t) correlacionadas espacialmente (Rouhani 
y Wackernagel, 1990; Rouhani et al., 1992, 
Knotters y Bierkens., 2001; Kyriakidis y Journel, 
2001a, 2001b). 

La metodología aquí propuesta corresponde 
al primer punto de vista, donde el proceso 
espacio-temporal se descompone en: 


Zlx,t)=m(x,t)+R(x,t) VxeD, VteT (3) 


en donde m(x,t) es una componente determi- 
nista conocida como “tendencia”. La variable 
aleatoria R(x,t), a la que llamaremos “residuo”, 
cumple con lo siguiente: 


a) Su esperanza matemática existe; no depen- 
de del punto de observación y es nula. 

b) Para todos los 
varianza finita. 


incrementos tiene una 


Como en el caso espacial, estas suposiciones 
se conocen como “hipótesis intrínseca” y es 
suficiente saber lo dos primeros momentos de 
la variable para caracterizarla. 


Sea Z(x,t) una VA Z(-) en la posición x=(x,y), 
al tiempo t. 

Entonces, el semivariograma muestral para 
los incrementos Ax=(Ax, Ay) y At se calcula con 
base en la siguiente expresión: 


y, (Ax, At) (4) 
O )-261)') 
= ZU + Ax, t¡+ At)-Z(x,t 
2(N(x t)) 2, 2, Ñ 
en donde __1__ son el número de pares 
2(N0s t)) 


separados por los incrementos Ax y At; por 
definición y(Ax=0, At=0)=0. 

Finalmente, la estimación de la carga 
hidráulica se realiza utilizando la técnica 
de krigeado espacio-tiempo (Ibañez, 2003; 
Mendoza, 2008). 


Metodología y resultados 


En esta sección se muestran los resultados. 
Por ello, consideramos conveniente hacer una 
descripción de la variable analizada desde los 
puntos de vista espacial, temporal y espacio- 
temporal. Se presentan sus estadígrafos con 
intención de que, al momento de modelar el 
semivariograma espacio-temporal, se tenga 
una idea de la forma en cómo varía en cada 
uno. 


Estadígrafos en espacio 


Como primer paso, seanalizó en formaseparada 
la información espacial de la temporal. De 
esta manera, se describieron los estadígrafos 
que caracterizan la carga hidráulica para cada 
año, utilizando las librerías desarrolladas por 
Deutsch y Journel (1998). 

En el diagrama de barras de la figura 2 se 
muestra el número de lecturas anuales de 
pozos de observación. Cabe destacar que se 
cuenta con información rica en el espacio, pero 
en algunos tiempos, con pocas observaciones. 
Para el periodo 1970-2004 existen años que van 
desde cero lecturas (2000) hasta un máximo de 
170 (2001). 
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La base de datos en ocasiones es rica en 
espacio, pero en otras muy pobre, lo que afecta 
directamente al cálculo de la estructura del 
semivariograma. Por tal motivo se procede a 
trabajar con aquellos años que presentan treinta 
o más observaciones, mismas que se muestran 
en la figura 3a. Se observa cómo la muestra se 
redujo a veinte columnas representantes de 
los años de muestreo, que cubren desde 1981 
hasta 2004, teniendo como número mínimo de 
lecturas treinta, para 1997, y máximo de 173, 
para 2001. 

Los estadígrafos que caracterizan las 
lecturas de la carga hidráulica se resumen 
en el cuadro 1. Los estadígrafos indican que 
en todos los años los valores de la media y 
mediana están separados, y se reportan valores 
de desviaciones estándar altas con valores de 
simetría que señalan sesgo en la muestra y 


forma platicúrtica. 


Por lo antes comentado, será necesario 
transformar los datos, con objeto de trabajar 
con muestras cuyos valores asemejen una 
distribución simétrica. 


Estadígrafos en tiempo 


A diferencia del análisis exploratorio espacial, 
el temporal toma en cuenta la variación de la 
variable aleatoria, carga hidráulica, en una 
misma posición; es decir, analiza la serie de 
tiempo de la carga hidráulica que se presenta 
en un pozo de observación. 

Observando la base de datos, es evidente 
que la información temporal no se encuentra 
densamente muestreada; basta considerar que 
de los 322 pozos existentes, para uno hay un 
máximo de 31 lecturas y, como mínimo, una. La 
figura 4 muestra el histograma del número de 
observaciones en el periodo 1970-2004. 
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Figura 2. Histograma número de pozos muestreados en el periodo 1970-2004. 
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Figura 3. Histograma, a) pozos con 30 observaciones; b) 13 pozos con 19 o más lecturas anuales. 
Cuadro 1. Estadígrafos de la carga hidráulica. 
Núm. obs. Media Mediana Varianza Desv. est. Mínimo Máximo Simetría Curtosis Fecha 
69.00 1767.91 1 752.68 1 998.79 44.71 1 690.95 1 913.70 1,91 3.38 01/08/1981 
67.00 1754.22 1740.13 2124.79 46.10 1 678.06 1 903.44 2.06 DN 01/08/1985 
83.00 1731.84 1720.70 2 498.33 49.98 1 654.49 1 897.99 1.82 3.53 01/08/1990 
38.00 1745.19 1716.31 4 680.44 68.41 1 642.32 1 920.48 1.36 0.90 01/05/1993 
85.00 1737.35 1 713.48 3377.93 58.12 1 640.42 1 894.43 1.60 1.64 01/12/1993 
33.00 1724.72 1711.25 2.338,95 48.36 1 635.99 1891.99 1.64 4.09 01/05/1994 
a 44.00 1 748.38 1713.01 5 109.31 71.48 1 632.49 1907.93 1.07 -0.07 01/12/1994 
a 49.00 1732.41 1 708.35 3 934.71 62.73 1 654.57 1907.73 1.73 1.95 01/05/1995 
= 112.00 1728.37 1706.89 2 965.89 54.46 1 629.24 1 892.48 1.84 2.72 01/12/1995 
E 33.00 1738.13 1705.61 4 669.56 68.33 1 656.41 1907.30 1.51 0.86 01/05/1996 
= 31.00 1 725.83 1 701.88 3 853.44 62.08 1 654.79 1 890.48 1.97 2.77 01/12/1996 
a 30.00 1 730.73 1705.16 3733.14 61.10 1 658.75 1 888.70 1.61 1.49 01/12/1997 
E, 40.00 1722.49 1700.56 2 554.13 50.54 1 655.65 1 888.08 2.02 4.15 01/12/1998 
ZÉ 75.00 1 720.65 1 693.65 3 862.62 62.15 1 623.33 1 890.40 ¿des 1.55 01/12/1999 
a 145.00 1704.16 1 686.73 4 554.91 67.49 1 466.44 1 910.58 1.14 2.49 03/12/2003 


« vj30]0u72], 


Mendoza-Cázares, EY. y Herrera-Zamarrón, G. del S., Estimación espacio-temporal de la carga hidráulica utilizando el concepto de... 


Por dicha razón, se decidió trabajar con 
pozos que reportaran, en el tiempo, 19 o más 
lecturas; así se trabajó con una submuestra que 
contabiliza un total de 13 pozos, tal como se 
observa en la figura 3b. Basta comentar que la 
decisión anterior obedece a que al calcular el 
semivariograma muestral en tiempo, serequiere 
de un número mínimo de observaciones 
para evitar que el semivariograma tenga un 
comportamiento errático. 

El cuadro 2 presenta los estadígrafos de los 
trece pozos. Hacemos notar que la muestra no 
guarda un comportamiento de una distribución 
simétrica. 


Estimación de la deriva 


Es común que los niveles piezométricos 
muestren una pendiente global en la dirección 
del flujo, por lo que suponemos que la media 
de la variable sí depende de su posición en el 
espacio y el tiempo; es decir, existe una deriva. 
Se obtuvieron los semivariogramas muestrales 
para los niveles piezométricos para los años 


con más datos y se observó que éstos crecen 
más rápidamente que h?. Se sabe que una 
función aleatoria intrínseca es aquella cuyos 
incrementos tienen esperanza matemática 
y varianza definidas e independientes de la 
posición para todo vector h (Samper y Carrera, 
1990), por lo que en este caso la variable es 
intrínseca y es lógico pensar que existe una 
derivada. 

Para quitar el efecto de deriva se procede 
a tratarla como indica la ecuación (3) y 
expresando la deriva como: 


mlxa,t)= Y bifelot) Vs eDxT 65) 
k=0 


donde m(x,,t) está compuesta por k funciones 
conocidas f(x,t) seleccionadas para ajustar 
el valor de deriva y b,, coeficientes del ajuste 
asociados con la k-ésima función. En este 
caso, el ajuste de la deriva fue a través de un 
polinomio de primer grado, como indica la 
ecuación (6): 
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Figura 4. Histograma número de observaciones por pozo en el periodo 1970-2004. 
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Cuadro 2. Estadígrafos de la lectura de la carga hidráulica por pozo. 
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Media Mediana Varianza Desviación Mínimo Máximo Simetría Curtosis Pozo 

1740.13 1747.87 721.01 26.85 1 698.59 1775.12 -0.183 -1.561 613 
17S8:21 1744.06 1 008.98 31.76 1 683.23 1 801.13 -0.050 -1.056 605 
1736.31 1740.77 570.89 23.89 1691.15 1769.08 -0.212 -1.300 612 
1 802.63 1798.99 163.58 1279 1780.13 1 827.29 0.397 -0.511 52 
1735.56 1734.91 535.23 23.14 1 689.01 1769.79 -0.158 -0.992 641 
1 743.29 1750.11 590.28 24.30 1 690.75 1770.8 -0.663 -0.822 0610-A 
1 891.61 1 892.85 135.64 11.65 1 858.03 1 910.92 -0.748 1.686 759 
1718.47 U7ALZALS) 268.92 16.40 1 700.88 1765.85 1.140 1.512 870 
1 900.33 1 899.14 86.55 9.30 1 886.82 1 912.29 0.078 -1.545 0548-A 
1470927 1710.74 529.00 23.00 1 677.69 M7 0.349 -0.830 629 
1741.17 1743.61 493.91 22.22 1 688.07 1769.17 -0.824 0.284 627 
1 700.39 1702.40 174.19 13.20 167985 1724.11 -0.108 -0.964 0999-A 
1705.10 1706.85 182.01 13.49 1679.12 1730.42 0.034 -0.722 1638 

mx, t)=a8,+0/x+0,y+0st (6) Análisis estructural 


donde los valores de los coeficientes son: 


a,=-4 268.5, a,=-0.000044, a,=0.002653, 
a,=-1.92210 


Una vez determinada la deriva, ésta se 
resta de la observación inicial y se trabaja con 
los residuos. Suponemos que los residuos son 
aleatorios y se comportan como funciones 
aleatorias intrínsecas. Los estadígrafos de 
los residuos, en el espacio, se muestran en el 
cuadro 3 y, en el tiempo, en el cuadro 4. 

En resumen, desde el punto de vista 
espacial, setienequela variable carga hidráulica 
presenta un comportamiento ascendente en 
la distancia, típico en acuíferos, que refleja 
un efecto de deriva. Se considera que los 
residuos son estacionarios de segundo orden 
(es decir, existe la media y no depende de la 
posición, y la covarianza existe y depende 
únicamente del vector de separación entre 
las posiciones), de esta forma se trabajará 
con ellos para determinar su estructura y, 
posteriormente, calcular los parámetros que 
los caracterizan. 


El objetivo de esta sección consiste en 
determinar la estructura de los datos en 
espacio-tiempo, por lo que se procederá a 
calcular el semivariograma muestral de los 
años identificados en el cuadro 3. Se calculó el 
semivariograma espacio-temporal y luego se 
efectuó la validación cruzada del modelo. 

Las propuestas de Bilonik (1983), Cameron 
y Hunter (2006), y Gardner y Sullivan (2004) 
se basan en calcular un semivariograma 
promedio, partiendo de cada semivariograma 
individual, en una distancia “lag” común. 
Este tipo de semivariograma es denominado 
“compuesto”. En este trabajo se propone 
analizar un semivariograma que de forma 
conjunta integre la información espacial y 
temporal, tal como se indica en la siguiente 
sección. 


Semivariograma espacio-temporal 


En esta sección se analiza la estructura del 
semivariograma espacio-tiempo, denominada 
“semivariograma de superficie”, y el ajuste de 
un modelo teórico producto suma al mismo. 
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Cuadro 3. Estadígrafos de los residuos de la carga hidráulica. 


Núm. obs. | Media | Mediana | Varianza | Desv.est. | Mínimo | Máximo | Simetría | Curtosis Fecha 
69.00 4.83 2.12 846.05 29.09 -70.22 78.20 0.22 1.14 Oct-81 
67.00 3.57 -0.61 987.37 31.42 -61.36 130.28 1.30 Sl Ago-85 
83.00 -5.18 -8.62 1162.79 34.10 -73.89 98.95 0.89 1.80 Ago-90 
38.00 -4.30 -4.09 1 028.02 32.06 -79.91 62.54 -0.06 0.38 May-93 
85.00 2.14 -2.58 1190.12 34.50 -80.52 106.37 0.83 1.57 Dic-93 
33.00 -6.79 -4.55 873.93 29.56 -84.02 62.71 -0.11 1.25 May-94 
44.00 -2.09 -5.48 1 432.86 37.85 -86.24 110.36 0.62 1.29 Dic-94 
49.00 -5.87 -6.49 989.02 31.45 -61.47 98.65 iL1l PRSil May-95 

112.00 1.83 -3.07 1071.28 32.73 -87.28 108.59 0.59 1.58 Dic-95 
33.00 Bi2Zo -4.64 1 088.46 S2I99 -59.10 103.34 1.68 3.82 May-96 
31.00 -8.34 -13.89 1321.18 36.35 -57.75 112.40 2.19 5.86 Dic-96 
30.00 -0.33 -11.39 1 931.63 43.95 -55.34 112.89 1.60 187 Dic-97 
55.00 0.31 -6.29 777.00 27.87 -56.04 78.89 0.68 1.01 Dic-98 
75.00 -1.46 -9.25 1 473.56 38.39 -84.33 118.24 0.99 1.63 Dic-99 

170.00 5.21 -1.67 1 324.63 36.40 -107.68 114.31 0.41 2.03 Dic-01 
44.00 2.64 3.36 1511.36 38.88 -109.71 90:09 -0.41 1.63 May-02 
44.00 6.92 8.83 1 462.39 38.24 -106.87 97.88 -0.50 1.99 Oct-02 
44.00 12.05 15.69 1 474.57 38.40 -104.98 87.38 -0.74 1.86 Nov-03 

145.00 -10.94 -14.89 2 270.32 47.65 -243.50 147.52 -0.11 1.49 Dic-03 
41.00 14.69 15.90 1 652.92 40.66 -105.14 11102 107 -0.26 11,55) May-04 

Cuadro 4. Estadígrafos de los residuos de la carga hidráulica por pozo. 
Media Mediana Varianza Desviación Mínimo Máximo Simetría Curtosis Pozo 

-31.38 -31.82 25.95 5.09 -37.94 -19.88 0.94 0.74 613 
-2.81 -2.36 81.39 9.02 -14.10 15.81 0.94 0.42 605 

2.81 4.04 35.57 5.96 -7.98 14.45 0.07 1.42 612 

59.84 62.66 -1.45 2.60 752 

-27.00 -24.53 -0.01 -1.51 641 
-9.24 -9.22 Lil) 2.08 0610-A 
-9.82 -8.95 -1.15 1.84 759 
-4.20 -4.32 1.92 5.62 870 

-21.20 -19.55 13.68 3.70 -27.89 -17.63 -1.10 -0.32 0548-A 

-11.79 -13.17 44.41 6.66 -20.10 3.86 0.97 1.10 629 

-13.88 -14.27 47.11 6.86 -23.37 -4.10 0.10 1.73 627 
-4.75 -5.44 20.73 4.55 -11.09 5.05 0.48 -0.11 0999-A 
-6.18 -4,47 9.67 3.11 -11.14 -2.55 -0.39 -1.76 1638 


Nota: debido a la decisión que se hizo en el espacio (tomar en cuenta los años donde se tienen treinta o más mediciones), el 
número de observaciones en el tiempo se ve afectada, ya que se eliminaron datos que van de 1971 a 1981, por lo que el número 
de observaciones para un mismo pozo se ven reducidas. Por el momento, los estadígrafos corresponden a los mismos pozos que 


se muestran en el cuadro 2. 
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Se procedió a calcular el semivariograma 
muestral (ecuación (4)), utilizando el algoritmo 
propuesto por De Cesare et al. (2002). El 
semivariograma espacio-temporal se presenta 
en la figura 5. 

De Cesare et al. (2001) proponen el modelo 
que combina el producto y suma de la siguiente 
manera: 


Cosllts, h,) > k¡C, (1,)C,(h1,) + k.C, (h,)+ kC; (1,) (7) 


o, de forma equivalente, el semivariograma: 


stas h,) => [x, +k,CC(0)ly, (1, ) 
+ [kz +k ¡E,Csto lv.) kv, (,) y, (1,) 


(8) 


donde C, y C, son la función de covarianza en 
espacio y tiempo, respectivamente; y, y y, son 
las correspondientes funciones de semivario- 
grama, en espacio y tiempo. Es de remarcar 
que C,(0) es la cima de y,; C(0) es la cima 
de y, y C/(0) es la cima de y, y por definición 
y, (0,0)=y.(0)=y,(0)=0. Por otro lado, se deben 
cumplir las condiciones siguientes: k,+k,C,(0)=1 
y k,+k,C.(0)=1. Para ajustar el modelo teórico al 
semivariograma muestral, se requiere estimar 
la cima global, espacio-tiempo, misma que se 
puede definir a través de un procedimiento 
de ajuste visual. En este estudio se evaluaron 
los modelos teóricos (esférico, exponencial, 
Gaussiano, potencia y agujero) presentes en las 
utilerías del GSLIB. 
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Figura 5. Semivariograma de superficie espacio-temporal, a) vista en la escala del espacio, 


b) vista desde la escala del tiempo, c) 3D. 
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El procedimiento para el ajuste del 
semivariograma espacio-temporal se puede 
hacer valiéndose de los semivariogramas 
teóricos en espacio y tiempo, y la selección de 
la cima global por medio de un ajuste visual 
que reproduzca el semivariograma muestral. 
Este proceso depende del analista, por lo que 
el valor obtenido varía dependiendo de quién 
lo proponga. Este proceso no es ajeno a la 
metodología que aplican los geoestadísticos en 
el caso de análisis univariados. 

Como 


primer paso, se 


para diferentes 


graficó el 
semivariograma muestral 
incrementos de distancia en el eje del espacio y 
en el eje para el tiempo. Efectuado lo anterior, se 
propuso un semivariograma, producto suma, 
donde se conoce la cima de la componente 
espacial y temporal. 

Para ello se propusieron diferentes valores 
de los parámetros presentes en la ecuación 
(8), hasta ajustar el semivariograma teórico 
al muestral, buscando que la media de las 
diferencias, que denominaremos “errores”, 
entre el valor del semivariograma teórico y el 
muestral, fuera mínima. Siguiendo los pasos 
anteriores, se llegó al semivariograma teórico 
mostrado en la figura 6: el valor de error es de 
12.59. Este proceso es de prueba y error hasta 
determinar el mejor ajuste; en el artículo no 
se mencionan los diferentes valores de error 


Semivariograma espacio-tiempo 


4 000 ] 


o Muestral 
a Teórico 


Figura 6. Semivariograma muestral y teórico 
espacio-temporal. 


ajustados para los diversos modelos teóricos. 

Al cuantificar el error de las mediciones 
entre la función de semivariograma y el 
respectivo modelo teórico, se obtiene una forma 
cuantitativa para seleccionar el mejor modelo 
que, en este caso, da como resultado un valor 
de la cima, C (0,0), igual a 1 840. 

Es importante remarcar que este criterio no 
es automático, ya que en el proceso de ajuste 
del semivariograma espacio-tiempo se van 
proponiendo diferentes valores hasta llegar a 
un valor que, de acuerdo con los errores, sea 
el menor. Sin embargo, dicho proceso permite 
al modelador ajustar de acuerdo con su criterio 
de mejor manera en las interdistancias de 
espacio O tiempo que muestren un mayor 
número de pares (ya que los valores de éstas 
en general debieran ser más representativas 
del semivariograma), y desechando las que 
no sean congruentes; esto se ve reflejado 
en la figura 6. El gráfico muestra un mejor 
ajuste en las interdistancias (h) más cercanas. 
Conforme se incrementa la interdistancia, el 
semivariograma muestral presenta mayor 
dispersión y el semivariograma teórico se ajusta 
en la porción central. Verificando el número 
de pares para los puntos con interdistancias 
mayores, se confirma que su valor es bajo. 
Este tipo de criterios son importantes para que 
el modelador, bajo su criterio y experiencia, 
decida el mejor ajuste. Esto no pasa con un 
proceso de ajuste automático, sea por mínimos 
cuadrados o mínimos cuadrados ponderados. 

Es importante destacar que los resultados 
se ven afectados al momento de decidir el valor 
de la cima del modelo espacio-tiempo. Isaaks 
y Srivastava (1989) ejemplifican la sensibilidad 
de los datos al seleccionar un modelo, tratando 
el efecto que tienen cuando se varía la escala 
(diferente cima), la forma (diferente modelos), 
pepita (diferente pepita), el rango (distintos 
rangos) y la anisotropía (distintas orienta- 
ciones). Cuando se varía la cima de los modelos, 
no existe un efecto en los valores estimados y 
los pesos del krigeado. Sin embargo, sí existe 
una afectación en la varianza, modificando el 
valor de la varianza en la misma proporción 
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(factor de escalamiento) en que se varía la cima 
del modelo propuesto. 

Los ejemplos de Issaks y Srivastava orientan 
de forma intuitiva al modelador acerca de 
la sensibilidad en la estimación de la carga 
hidráulica y el valor de la varianza, al variar 
los parámetros del modelo. 

Por esta razón es necesario llevar a cabo un 
análisis comparativo. En este trabajo se hizo a 
partir del análisis de los errores de la validación 
cruzada, lo cual se comenta en la siguiente 
sección. 


Validación cruzada 


Para verificar qué tan bien 


semivariograma teórico, se procede a utilizar 


ajusta el 


el método de validación cruzada para después 
calcular el error medio (EM), el error cuadrático 
medio (ECM) y error cuadrático medio estándar 
(ECME) (ver la expresión matemática de cada 
uno de ellos en el cuadro 5). 


El EM mide la 
subestimación de la variable; si es positivo, 


sobreestimación O 


la variable en promedio se subestima y si es 
negativo lo contrario; la estimación es mejor si 
éste es cercano a cero. El ECM da una medida 
promedio de la magnitud al cuadrado de los 
errores; también se espera que el valor sea 
pequeño. Con respecto al ECME, los valores 
se interpretan como el ajuste que presenta 
la varianza del modelo del semivariograma, 
espacial y espacio-temporal, con los errores 
de estimación; por esta razón, cuando el valor 
del ECME es cercano a la unidad, se interpreta 
como que la varianza del modelo teórico y el 
error cuadrático de la estimación en promedio 
son cercanos. 

El proceso de validación cruzada se realiza 
utilizando el algoritmo propuesto por De 
Cesare et al. (2002). Con el objetivo de evaluar 
los resultados, se seleccionaron los años 1993, 
1995 y 1999, que corresponden a los mismos 
años estimados con el método multivariado de 


Cuadro 5. Distribución de errores utilizando el método de validación cruzada con krigeado ordinario (KO), con estimación 


espacio-temporal cokrigeado (Cok) y espacio-temporal de variables aleatorias (ET). Mín = mínimo; máx = máximo; 4 = media; 


n n 
0 varianza; EM = error medio = 1 i2(x)-Z 5) ; ECM = error cuadrático medio = - NilzG)-z JD y ECME = error 
¡=1 


$ 2 
1 q ' 
cuadrático medio estándar = a ge = (1) 


i=l 


n 
i=l o; 


Valores estimados de los residuos, con el método de a 

Método y año validación cruzada 
Mín Máx p o? EM ECM ECME 
KO, 1993 -60.192 60.829 -0.192 418.360 -0.190 390.550 0.820 
KO, 1995 -66.762 72.764 -0.098 463.750 -0.100 252.340 0.800 
KO, 1999 -57.691 59.654 -0.396 446.104 -0.400 426.300 0.800 
Cok., 1993 -89.733 69.223 -0.815 599.513 -0.790 349.510 0.920 
Cok., 1995 -89.081 78.396 -0.620 611.010 -1.630 284.310 0.990 
Cok., 1999 -85.555 71.300 0.071 723.890 -1.960 309.030 0.950 
EIOOS -65.868 98.995 1.018 998.881 -1.045 95595 0.408 
ET, 1995 -57.191 100.489 1.699 797.711 -0.131 137.643 0.361 
EOS -54.158 96.339 2.426 877.105 -0.967 339.706 0.683 
ET 1981-2004 -131.265 122.227 -2.014 878.985 -0.621 318.731 0.804 
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funciones aleatorias (Mendoza y Herrera, 2007). 
La comparación entre los resultados obtenidos 
para los métodos de funciones aleatorias y 
multivariado se discute más adelante. 

Los resultados de validación cruzada para 
1993 muestran que los valores de error son los 
siguientes: EM=-1.045, ECM=195.53, ECME= 
0.408; el valor promedio de desviación estándar 
es de 31.60. La figura 7a muestra el diagrama 
de dispersión entre los valores observados de 
la carga hidráulica para 1993 y los estimados 
por el método de validación cruzada. 

Los resultados de validación cruzada 
para 1995 muestran que los valores de error 
son los siguientes: EM=-0.131, ECM=137.64, 
ECME=0.361; el valor promedio de desviación 
estándar es de 18.0. La figura 7b muestra 
el diagrama de dispersión entre los valores 


observados de la carga hidráulica para 1995 
y los estimados por el método de validación 
cruzada. 

Por último, los resultados de validación 
cruzada ¡para 1999 muestran siguientes 
valores de error: EM=-0.967, ECM=339.706, 
ECME=0.683; el valor promedio de desviación 
estándar es de 29.61. La figura 7c muestra 
el diagrama de dispersión entre los valores 
observados de la carga hidráulica para 1999 
y los estimados por el método de validación 
cruzada. 


Estimación espacio-temporal 


Los resultados de la estimación de la carga 
hidráulica en el acuífero de Querétaro- 
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Figura 7. Diagrama de dispersión; a) carga hidráulica del año 1993, b) carga hidráulica del año 1995, 
c) carga hidráulica del año 1999. 
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columnas y 51 renglones, con una distancia 
entre renglón y columna de 1 100 m y un delta 
de tiempo de una unidad. Sólo se reportan 
los planos de la carga hidráulica de los años 
1993, 1995 y 1999, dejando de lado los años 
intermedios y la predicción correspondiente 
a los años posteriores a 2003, último año que 
cuenta con lecturas. 

Las ilustraciones 8 a 13 muestran los mapas 
dela carga hidráulica y la varianza paralos años 
comentados, y de forma global se describen los 
siguientes patrones: en la figura 8 se muestra 
la estimación de la carga hidráulica para el año 
1993; se observa que las líneas equipotenciales 
son acordes con las mediciones reportadas 
para cada pozo de observación. Cabe destacar 
que se marca un parteaguas en la zona de la 
laguna El Salitre y en el área del acuífero de 
Querétaro se tiene una mejor definición de las 
zonas con mayor abatimiento. Esta definición 
se observa en el mapa como áreas concéntricas 


de una influencia local que difieren del entorno. 
Al verificar estos puntos corroboramos que, 
efectivamente, corresponden a pozos donde 
existe una carga hidráulica menor al entorno. 

En el caso de La Cañada y zonas altas de 
la cuenca, se observa el valor subestimado 
de la carga hidráulica; sin embargo, de forma 
general, se mantiene el mismo patrón de las 
líneas de flujo en dirección al centro del valle 
de Querétaro. 

El plano de varianza (figura 9) presenta que 
la porción central del acuífero está mayormente 
muestreada; la línea equipotencial de 400 
unidades se distribuyen en un amplio sector y, 
las márgenes del acuífero pierden esta densidad 
de muestreo. 

Respecto a la estimación proyectada en el 
plano para 1995 (figura 10), se mantiene quelas 
estimaciones son acordes con las mediciones 
reportadas para cada pozo de observación. 
El parteaguas en la zona de la laguna El 
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Figura 8. Krigeado espacio-temporal para el año 1993, estimación de la carga hidráulica. 
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Figura 9. Krigeado espacio-temporal para el año de 1993, valores de varianza. 
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Figura 10. Krigeado espacio-temporal para el año 1995, estimación de la carga hidráulica. 
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Salitre se mantiene, así como la definición de 
las zonas con mayor abatimiento. En el caso 
de La Cañada y zonas altas de la cuenca, se 
observa que el valor de la carga hidráulica es 
subestimado; sin embargo, de forma general 
se continúa el mismo patrón de las líneas 
de flujo en dirección al centro del valle de 
Querétaro. 

En la figura 11 se muestra el plano de 
varianza, que mantiene el mismo patrón del 
año 1993, marcando una diferencia en la zona 
de la laguna El Salitre, producto del incremento 
de pozos de muestreo para el año 1995. 

La estimación proyectada en el plano para 
1999 (figura 12) presenta una correspondencia 
entre lo estimado y las mediciones reportadas 
para cada pozo. Cabe señalar que para este 
año no se marca un parteaguas en la zona 
de la laguna El Salitre, figurando que los 
acuíferos Los Apaseos, Obrajuelo y Querétaro 


han sido afectados en su dinámica original 
por su inversión en el flujo, sobre todo en 
Los Apaseos, ahora centro del acuífero del 
valle de Querétaro. Esta hipótesis deberá 
ser corroborada con información adicional 
obtenida para esta zona. Este parteaguas 
es acorde con el modelo conceptual 
comentado en el área de estudio para el 
caso no estacionario, pero la posición del 
mismo no está bien definida y mucho menos 
confirmada. 

Que el valor del error cuadrático medio 
estándar sea elevado; de hecho, el más elevado 
de los tres años analizados, indica de forma 
tácita que el modelo del semivariograma 
espacio-tiempo no es bueno para estimar este 
año en particular. Es de esperar que en los 
puntos donde se estima, dentro de la malla 
indicada, no sean siempre acordes con la 
realidad. 
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Figura 11. Krigeado espacio-temporal para el año de 1995, valores de varianza. 
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Figura 12. Krigeado espacio-temporal para el año 1999, estimación de la carga hidráulica. 
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Figura 13. Krigeado espacio-temporal para el año de 1999, valores de varianza. 
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El plano de varianza (figura 13) presenta 
un aumento de la varianza en la zona oeste, 
producto de una reducción de los pozos de 
muestreo. 

Estimación en tiempo 


En este punto se muestran los resultados 
obtenidos con el método de validación 
cruzada en los pozos que cuentan con 19 o más 
observaciones de la carga hidráulica, mismos 
que se listan en el cuadro 4. La intención 
es mostrar las lecturas de estimación de la 
carga hidráulica, comparadas con los valores 
observados, cuando se ha quitado el dato que 
se estima. De esta forma, sólo se mostrarán los 
resultados para luego ser discutidos. La figura 
14 muestra tres hidrógrafos de los pozos 759, 
752 y 629. 

La espacial de los pozos 
seleccionados se muestra en la figura 15. Es fácil 


posición 


identificar que los aprovechamientos no están 
en una zona restringida; por el contrario, se 
ubican en todo el acuífero. De forma particular, 
el pozo 752 se sitúa en la zona de recarga; el 
pozo 629, en la porción central del acuífero 
Querétaro y el pozo 759, en la parte alta del 
acuífero. La falta de información histórica en 
la porción del acuífero Obrajuelo se da porque 
el mínimo de observaciones, 19 lecturas desde 
el año 1981 hasta el 2004, no se cumplió. Esta 
imposición obedece a que es necesario contar 
con un número mínimo de observaciones para 
generar los semivariogramas en el tiempo. 
Aclarado el punto anterior, se comentan los 
resultados de los hidrógrafos. El pozo 759 se 
ubica en la porción alta del acuífero y de manera 
general describe un descenso progresivo de la 
carga hidráulica. Para el año de 1981 se reporta 
una elevación de la carga hidráulica de 1 910.92 
m y, para el año del 2003, se reportan 1 879.14 


Elevación de la carga hidráulica, pozo 759 


Elevación de la carga hidráulica, pozo 752 
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Figura 14. Hidrógrafos pozos 759, 752, 629. 
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Figura 15. Pozos de observación con más de 19 lecturas, en el tiempo, de la carga hidráulica. 


m; la diferencia en estos 23 años es de 31.78 
metros. 

El descenso de la carga hidráulica es 
constante, de forma gradual, existiendo 
algunos puntos de recuperación, como se 
distingue en el año de 1990, donde el nivel se 
recupera 3.0 m, tomando en cuenta la tendencia 
general, para después descender de manera 
gradual. Cabe resaltar que en el tiempo, entre 
los años 1999 y 2001, se mantiene constante, 
para después descender 6.51 metros. 

Respecto a los valores estimados, se 
observa que de manera general se ajustan a la 
tendencia, presentando las mayores diferencias 
en las series donde no existen observaciones 
en el año inmediato anterior y posterior (ver 
periodo 1981-1990), pero en los periodos donde 
hay observaciones, las estimaciones presentan 
un mejor ajuste; de hecho, las diferencias entre 


ambas no sobrepasan los 1.5 metros. 


El pozo 752 se ubica en la zona de recarga 
del acuífero. De manera general, describe un 
descenso progresivo de la carga hidráulica. 
Para 1988 se reporta una elevación de la carga 
hidráulica de 1 796.45 m y, para el año 2003, 
se reportan 1 788.43 m. La diferencia en estos 
16 años es de 8.02 m. El descenso de la carga 
hidráulica es gradual. 

Respecto a los valores estimados, es 
evidente que la estimación de la carga 
hidráulica no se ajusta a las mediciones en 
los pozos que se encuentran en las márgenes 
del acuífero. En estos sitios se identifica que 
las estimaciones están por debajo de las 
mediciones; es decir, se subestiman. La falta 
de información, en el tiempo, ocasiona que 
los puntos de estimación sean extrapolados. 
Los valores de varianza indican que las 
estimaciones presentan una incertidumbre 
alta. 
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Los valores estimados, que presentan una 
mayor diferencia entre los valores medidos, se 
identifican en los intervalos de tiempo donde 
no existen datos antes o después. Lo contrario 
sucede en los intervalos de tiempo donde hay 
información en tiempos cortos. 

Journel y Rossi (1989) discuten el efecto 
que tienen las estimaciones, en los casos 
dónde se extrapola; aunque es aplicado en 
el caso espacial, el efecto es el mismo en el 
tiempo. 

El pozo 629 se ubica en la porción central del 
acuífero Querétaro. De forma general describe 
un descenso progresivo de la carga hidráulica. 
Para 1981 se reporta una elevación de la carga 
hidráulica de 1 753.12 m, y para el 2003 se 
reportan 1 677.69 m. La diferencia en estos 23 
años es de 75.43 m. El descenso de la carga 
hidráulica es gradual. Respecto a los valores 
estimados, se observa que de manera global 
se ajustan a la tendencia, que es lo buscado 
cuando no se cuenta con observaciones. 


Discusión 


Recapitulando, tenemos que en el aspecto 
espacial se seleccionaron 23 años (cuadro 1), los 
cuales cumplen con treinta o más observaciones 
(pozos) en cada muestreo. Este criterio se 
contempla considerando que al determinar 
la estructura de la variable se requiere de 
un número mínimo para efectuarlo. En la 
literatura, para calcular el semivariograma, 
se recomienda utilizar un número mayor de 
puntos de observación (cien observaciones); 
sin embargo, en la práctica, se trabaja con lo 
que se tiene, siempre manteniendo en mente 
que las decisiones que se toman tienen un costo. 
En el aspecto temporal, se seleccionaron los 
pozos que presentan 19 o más observaciones, 
dando como resultado 13 pozos. 

Así, se procedió a efectuar el análisis 
exploratorio de los datos, el análisis estructural 
espacio-tiempo y la validación del modelo 
propuesto, por lo que los comentarios de la 
discusión se enfocan hacia la validación y la 
estimación. 


Es importante destacar que la validación 
cruzada se llevó a cabo tomando en cuenta 
el modelo espacio-temporal, por lo que los 
resultados comprenden del año 1981 al 2004. 
De este periodo se seleccionaron los años 1993, 
1995 y 1999 para mostrar el grado de ajuste que 
tienen las estimaciones con las observaciones. 
El lector debe considerar que la estimación con 
este método no sólo muestra resultados en el 
espacio, sino también en el tiempo. 

Se optó por seleccionar estos años debido 
a que cuentan con un número importante de 
lecturas: 85, 112 y 75, respectivamente, lo cual 
favorece la estimación del semivariograma 
y, asimismo, para efectos de comparar los 
resultados con la técnica multivariada de 
funciones aleatorias. 

La figura 7a muestra el diagrama de 
dispersión entre los valores medidos de la 
carga hidráulica y los estimados para 1993. Los 
valores de la carga hidráulica medidos versus 
los estimados presentan un índice de regresión 
de R? = 0.94, que indica que éstos presentan 
una buena correlación. 

El cuadro 5 presenta un resumen de los 
errores de la validación cruzada obtenidos 
con los métodos kriging ordinario, cokriging 
y estimación espacio-tiempo. Una descripción 
detallada del procedimiento utilizado en los 
dos primeros métodos se puede encontrar 
en Mendoza y Herrera (2007). Para el caso 
espacio-tiempo de las variables aleatorias 
se utiliza la anotación ET para distinguirla 
de la técnica multivariada (Cok), también es 
importante explicar que cuando en el cuadro 
el método ET refiere a un año específico, se 
reportan los resultados de este método para 
el tiempo señalado (1993, 1995, 1999); es decir, 
la técnica estima en el espacio y también en el 
tiempo, utilizando los datos de todo el periodo, 
pero los estadígrafos mostrados en el cuadro se 
calcularon usando únicamente los datos para el 
año especificado. En contraste, lo que se reporta 
en la cuarta línea, con la anotación ET 19981- 
2004, son los estadígrafos de la validación 
cruzada manejando toda la información; o sea, 
tomando en cuenta los errores de la validación 
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de la carga hidráulica para el periodo 1981- 
2004. 

Aclarado lo anterior y retomando los 
resultados que muestra el cuadro, se identifica 
que los valores estimados de los residuos con 
el método de validación cruzada, para todos 
los casos, presentan un comportamiento 
simétrico. 

Analizando los valores del error (EM) para 
el caso de KO, se encuentra que son reducidos, 
cercanos a cero y todos ellos negativos, lo que 
indica una subestimación de los valores de la 
carga hidráulica. 

En el caso del Cok, la magnitud del EM es 
mayor que la del KO y, en todos los casos, es 
negativo, de la misma forma que para el KO. 

Para el ET, sólo en el caso donde se 
consideran todas las observaciones, el valor del 
EM es bajo, pero para cada año en particular, el 
error se incrementa un poco. 

En términos generales, el EM indica que 
para cada técnica, las estimaciones están siendo 
subestimadas. 

Para el caso del ECM, pensando éste como 
un valor de dispersión, ya que es el cuadrado 
de los errores, se observa que tanto en KO y 
Cok los valores son similares; sin embargo, 
al compararlos con el ET, el ECM se reduce 
drásticamente para los años 1993 y 1995; 
en 1999 se mantiene. El promedio del error 
cuadrático medio para los tres años en el ET 
es de 224.29 m, mientras que en el KO es de 
356.4 y en el Cok de 314.27. Esta reducción en 
el promedio del ECM para el ET se debe a que 
la estimación se realiza tomando en cuenta 
un número mayor de puntos de observación, 
no sólo en el espacio, sino también en el 
tiempo. 

De esta forma, el ECME indica que el mejor 
ajuste entre el error cuadrático y la varianza 
lo ofrece el método de Cok, seguido del KO, 
el ET para el periodo 1981-2004 y, por último, 
el ET para cada año. El hecho de que el ECME 
del ET para todo el periodo tenga un valor de 
0.804, mucho más cercano que los valores para 
cada año por separado, indica que la varianza 
del modelo teórico representa de mejor forma 


el error cuadrático para todo el periodo que 
para los años por separado. 

Es importante hacer notar que el ECM del 
ET indica que este método proporciona la 
representación de la carga hidráulica que mejor 
ajusta los valores observados; esto obedece a 
que el ET utiliza mayor número de puntos en la 
estimación, de igual forma que el Cok, sólo que 
el ET estima en los valores espaciales que se 
tienen en un año dado y utiliza la información 
temporal de todo el periodo; por esta razón el 
ECM es menor. 

Respecto a la estimación espacio-temporal 
proyectada en el plano, dejando al tiempo 
invariante, se observa que las estimaciones 
son acordes con las mediciones reportadas 
para cada pozo de observación. Cabe destacar 
que el método distingue un parteaguas en la 
zona de la laguna El Salitre y en el área del 
acuífero de Querétaro, y define las zonas con 
mayor abatimiento. Esta definición se observa 
en los planos como áreas concéntricas de una 
influencia local que difieren del entorno. 
Al verificar estos puntos corroboramos 
que, efectivamente, corresponden a pozos 
donde existe una carga hidráulica menor al 
entorno. 

En el caso de La Cañada y zonas altas de 
la cuenca, se observa el valor subestimado 
de la carga hidráulica; sin embargo, de forma 
general, se mantiene el mismo patrón de las 
líneas de flujo en dirección al centro del valle 
de Querétaro. 

Ahora bien, verificando los puntos que 
de acuerdo con los diagramas de dispersión 
(figura 7a) se alejan de la línea de 45 grados, 
correspondientes al caso donde la estimación es 
igual a la medición, tenemos que existen cinco 
valores: 1 721.94; 1 734.17; 1 766.72; 1 866.80 
y 1 890.94, que se ubican en las márgenes del 
acuífero (figura 8), caracterizados por pocos 
datos al hacer la estimación. Por la posición, 
se identifica que corresponden a puntos 
extrapolados. Verificando el valor de varianza 
que se obtiene para cada uno de estos puntos 
(figura 9) es posible corroborar que la varianza 
es alta. 


Tecnología y 


Ciencias del Agua, abril-junio de 2010 


Ciencias del Agua, abril-junio de 2010 


(vj30]0u32], 


Mendoza-Cázares, EY. y Herrera-Zamarrón, G. del S., Estimación espacio-temporal de la carga hidráulica utilizando el concepto de... 


En el caso de la piezometría para 
1995, se tiene que los valores de la carga 
hidráulica medidos versus los estimados, 
presentan un índice de regresión de 
R2=0.95 ( figura 7b), considerándose como 
buena correlación. 

Verificando los puntos que se resaltan 
en el diagrama de dispersión (figura 7b), 
tenemos que existen siete valores: 1 629.24; 1 
658.79; 1 715.60; 1734.20; 1 866.98; 1 874.85 y 1 
878.47, que al colocarse en el plano se ubican 
en las márgenes del acuífero (figura 10) y se 
caracterizan por tener pocos datos cercanos. 
Por la posición, al realizar la estimación, se 
identifica que no corresponde a un punto de 
interpolación; por el contrario, es un punto 
extrapolado. 

Verificando el valor de varianza que se 
obtiene para cada punto (figura 11) es posible 
corroborar que la varianza es alta. 

Al verificar los puntos que resaltan en el 
diagrama de dispersión (figura 7c) existen 
siete valores: 1 623.33; 1 650.5; 1669.5;1702.1; 
1720.05; 1792.06 y 1 864.61, que al colocarse 
en el plano se ubican en las márgenes del 
acuífero (figura 12), caracterizados por 
tener pocos datos para la estimación. 
Por la posición, al hacer la estimación, se 
identifica que no corresponde a un punto de 
interpolación; por el contrario, es un punto 
extrapolado. 

Verificando el valor de varianza que se 
obtiene para cada punto (figura 13) es posible 
corroborar que la varianza es alta. 

Con respecto a la varianza, es necesario 
resaltar que para todos los casos los planos son 
menores que los obtenidos con el método de 
krigeado. Esto es congruente con la teoría, ya 
que al aumentar el número de observaciones 
se reduce la varianza de la estimación. 

En los planos de varianza, se observa 
que en el acuífero Querétaro-Obrajuelo ésta 
es del orden de trescientas unidades que, al 
compararse con la obtenida con krigeado 
ordinario, se tiene que ha reducido en 
trescientas unidades, así como en la extensión 
del área de influencia. 


Conclusiones 


La estimación espacio-tiempo es una 
herramienta útil que puede ser tomada como 
alternativa para efectuar estimaciones de la 
carga hidráulica en un contexto espacial y 
temporal. 

La aplicación de esta técnica para 
determinar la carga hidráulica es innovadora, 
ya que resuelve el problema de los modelos 
lineales. La metodología desarrollada ofrece 
una alternativa para tratar las variables que 
presentan una deriva espacio-tiempo, así como 
indicar la forma de ajustar el modelo producto 
suma, propuesto por De laco et al. (op.cit.). 

Para realizar la verificación del modelo, se 
utilizó el método de validación cruzada que 
permite evaluar la certidumbre del modelo 
propuesto. Los valores de EM, ECM y ECME 
indican que el modelo producto suma se ajusta 
de manera adecuada a la información espacio- 
tiempo; sin embargo, al evaluar de forma 
individual los resultados de validación cruzada 
para los años seleccionados, se reportan algunos 
puntos con errores grandes. 

Las estimaciones de la carga hidráulica para 
los años 1993, 1995 y 1999 son congruentes con 
los valores observados en campo, así como la 
definición del parteaguas entre los acuíferos. 

Una ventaja del método espacio-tiempo es 
el incremento de observaciones para realizar 
la estimación; esta característica permite hacer 
estimaciones en zonas donde no se cuenta con 
información. 

El método es poderoso, ya que también 
realiza estimaciones en diversos tiempos. 
Esto se puede apreciar en la figura 14, donde 
se Observa que el método permite efectuar 
interpolaciones entre años donde se carece de 
información. 

El método espacio-tiempo requiere de 
un análisis estructural adecuado, ya que a 
partir de éste se seleccionan los parámetros 
que lo caracterizan. En nuestro caso se aplica 
una metodología en donde se contempla la 
selección de los parámetros, tomando en cuenta 
diferentes secciones de tiempo y espacio, y se 
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refuerzan con el índice de error. Esta propuesta 
es diferente a la de De Cesare et al. (op. cit.), ya 
que de manera gráfica se determina la cima 
del semivariograma de superficie. Así, ellos 
proponen un modelo producto suma que 
garantiza que el modelo es positivo definido y 
que no afecta en el sistema de krigeado. Estas 
dos innovaciones superan la propuesta inicial 
de Rouhani y Hall (1989). 

Finalmente, en términos de varianza, el 
método espacio-tiempo reduce el valor del 
error en la estimación; esto se puede observar 
comparando los resultados de varianza 
obtenidos con los métodos de krigeado 
univariado y krigeado multivariado; sin 
embargo, aun aumentando el número de 
puntos para hacer la estimación, se tienen 
zonas donde ésta no corresponde con las 
observaciones, ya que se ubican en la periferia 
del acuífero. Lo anterior obedece a que en la 
periferia las estimaciones son producto de la 
extrapolación que, aunado a los resultados del 
validación cruzada, son las áreas con mayor 
incertidumbre. Esta es una ventaja más de 
este tipo de métodos, ya que proporcionan 
información con la cual se pueden tomar 
confiabilidad 


de la estimación y las zonas de mayor 


decisiones considerando la 


incertidumbre. 
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Abstract 


MENDOZA-CÁZARES, E.Y. é HERRERA-ZAMARRÓN, G.S. Spatiotemporal estimation 
of hydraulic head using a single spatiotemporal random function model. Water Technology 
and Sciences, formerly Hydraulic engineering en Mexico (in Spanish). Vol. I, no. 2, 
April-June, 2010, pp. 87-111. 


This paper presents a geostatistical method for performing estimates using space-time kriging 
which is applied for the estimation of the water level of the Queretaro-Obrajuelo aquifer in the 
1981-2004 period. The estimates obtained by this method for the years 1993, 1995 and 1999 
are compared with ordinary kriging and cokriging using cross validation. The average mean 
error (ME) for the three chosen years was lowest for ordinary kriging (-0.23), the average 
mean squared error (MSE) is the lowest in the case of the space-time method (224.29) and the 
value of the average squared mean standard error (S5MSE) is better for the method of cokriging 
(0.95), because the SMSE values are close to unity. The SMSE for the space-time method is 
0.8 when considering all the times, but it decreases in particular for the selected years. The 
results of the estimated variances are always smaller using the space-time method, because it 
uses more information to get the estimate. Also, it is possible to make estimates in the whole 
space for all times. Therefore, it was concluded that the tool is powerful, since it considers all 
available information to produce the estimates. 


Keywords: single spatio-temporal model, geostatistics, hydraulic head estimation. 
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